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Non-perturbative renormalization in lattice QCD Yasumichi Aoki 

1. Introduction 

In many applications of lattice QCD it is necessary to deal with renormalization. A simplest 
example is to obtain the quark masses which are the fundamental parameters in the standard model. 
A bare quark mass is logarithmically divergent in the continuum limit, hence, is not defined without 
proper regularization and renormalization. A less simple example is to calculate a transition or 
decay amplitude of hadrons induced by interactions in the electroweak theory. The low energy 
effective operators are constructed through operator product expansion. They, in turn, are plugged 
in the QCD Lagrangian, with the Wilson coefficients calculated through (continuum) perturbation 
theory. The coefficients are often divergent, thus should be defined through renormalization with a 
certain renormalization scheme and scale. The lattice part is to calculate the matrix element of the 
operator between hadron states. The matrix element of the bare operator has the same divergence 
as the Wilson coefficient. In divergent cases, the matrix element or the operator must undergo 
the renormalization on the lattice. The same renormalization scheme and scale as the Wilson 
coefficient must be used to compensate the dependence on them, so that the final, physical matrix 
element is scheme and scale independent. 

As these typical examples show, there are two necessary steps required for the lattice renormal- 
ization: 1) removing the ultraviolet divergence in the observables and 2) matching to the scheme 
convenient to the continuum renormalization, such as MS. The step 2) can be performed in con- 
tinuum perturbation theory (cPT) once 1) is defined in a regularization independent (RI) way. The 
lattice perturbation theory (LPT) can be used to perform 1). However, the convergence of the 
series in general is not good due to tadpole contributions. Although a cure is provided by mean 
field improved perturbation theory, one has to deal with the ambiguity due to the choice of mean 
field coupling. In addition, for some quantities the matching 2) has been calculated to three loops 
in cPT. Such a calculation cannot be fully exploited if the LPT is used for the step 1) due to the 
difficulty of calculating two or more loops in LPT. These problems can be bypassed by the use of 
non-perturbative renormalization (NPR) applied for the step 1). 

The NPR on the lattice was born early in the previous decade. There have been dedicated 
plenary presentations in this series of conferences [|I], |[, Isj, ^. 

The popular NPR schemes to date are the momentum subtraction (RI/MOM) and Schrodinger 
functional (SF) schemes. The SF scheme often refers to its combined use with step scaling. The 
RI/MOM scheme is widely used in the light quark sector of state-of-the-art lattice calculations 
today. The use of SF scheme for light quarks is slightly restricted, presumably due to the need of 
tailored gauge configurations and subtleties when applying to the chiral fermions. The SF scheme 
has been extended to incorporate the heavy quark effective theory (HQET) on the lattice. It has 
been discussed in [Q] and more recently to some detail by Delia Morte [|5|]. 

As is well known there is the so called "window problem" in the RI/MOM scheme. To have 
minimal effect from the problem some practical prescription is required for the RI/MOM scheme. 
As an extension of the RI/MOM scheme, the RI/SMOM scheme was proposed recently and ap- 
plied to some realistic calculations. The new scheme is meant to ease the "window problem" by 
squeezing out unwanted non-perturbative effects from the setup of the scheme. This review focuses 
mainly on the recent development in RI/MOM and RI/SMOM schemes on how to work around the 
problem. 
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In what follows, RI/MOM and RI/SMOM schemes are discussed in Sec. 2. Recent develop- 
ments in SF scheme are given in Sec. 3. Sec. 4 contains application specific developments, such as 
the calculation of the strong coupling constant and quark masses. Concluding remarks are given in 
Sec. 5. 

2. Developments in the RI/MOM scheme 

The RI/MOM scheme [^ 0] is a popular non-perturbative renormalization scheme for multi- 
quark operators today. It has been applied to virtually all types of light quark discretizations used 
to date. The renormalization condition is imposed on the forward vertex function of the operator 
with external off-shell quark states having momentum p. This scheme is handy as no special 
gauge configurations are required. A drawback is that there is the so called window problem, 
i.e. the renormalization scale jx = ^J~p^ needs to be in the range AgcD ^ jU ^ 1/a. Finer lattice 
simulations are less problematic. However for the typical lattice cutoff used today, this could be a 
source of a sizable systematic error. In principle, the RI/MOM scheme can be used to determine the 
coefficient of the higher dimensional operators that appear in the Symanzik on-shell improvement. 
However, more operators have to be introduced than the SF scheme to properly disentangle the 
linear relation due to the use of off-shell external states, which would make the tuning difficult in 
practice. Operators with chiral or twisted mass fermions have no problem with this point since 
there are no 0{a) operators to mix at on-shell for typical applications. 

Several remarkable developments in the RI/MOM scheme and its extension are reported by 
this year. To discuss them it is required to have a review of the RI/MOM scheme. 

2.1 RI/MOM formulation 

For illustrative purpose let us go through the RI/MOM renormalization of a flavor non-singlet 
quark bilinear operator Or = uTd. The concept is easily extended to three, four • • • quark operators. 
Later on, the pseudoscalar/scalar operator is examined in detail for the quark mass renormalization. 
The bare operator Or with the Dirac structure F is renormalized with Zr as 

Of = ZrOr. (2.1) 

The RI/MOM renormalization condition is imposed to the amputated vertex function Tlr{p\,P2) 
with Landau-gauge fixed external mass-less off-shell quark states with the same incoming and 
out-going momenta pi = p2= p, 

^Tr[nr-^r] = l, (2.2) 

where Zq is the quark wavefunction renormalization factor, i.e. S^ = ZqS with S^ and S being the 
renormalized and bare quark propagator. The projection operator ^r is chosen so that the condition 
properly removes all the ultraviolet divergence and holds trivially at tree level with Zy = Zq = 1. 
The renormalization scale is set through the external off-shell momentum jx^ = p^. 

There are two different schemes due to slightly different definitions of the wavefunction renor- 
malization. The RI/MOM Zq is defined through the vector vertex with .'^y °^ 7^ 

|>|nv,.r„l = i. (2.3) 
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In the continuum or by the use of conserved vector current on the lattice, since Zy = 1 this deter- 
mines Zq. With the local vector current on the lattice Zy needs to be determine by other means. 
In the continuum, this condition is equivalent through the Ward-Takahashi (WT) identity with the 
condition on the inverse quark propagator. 



1 1 

Zqn 



-Tr 
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1. 



(2.4) 



Another scheme used to date is RI'/MOM scheme where instead of Eq. ( |2.3[ ) or ( [2.4| ) Zq is defined 
through 
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p-=^^ 



These two definitions of Zq are quite similar, which can be seen by the fact that the difference only 
appears at two loops in continuum perturbation theory [||]. To date RI'/MOM scheme is being used 
mostly for Wilson type fermions, while the RI/MOM is being utilized in domain-wall (DW) and 
overlap fermions. 

Once the Zy in RI/MOM or RI'/MOM scheme is obtained, it can be converted to that in other 
schemes by multiplying a conversion factor calculated in perturbation theory. 

These are the essence of RI/MOM schemes. In the lattice computation, there emerges a prac- 
tical issues associated with the use of momentum for the renormalization scale. The so-called 
window problem arises as the renormalization scale needs to satisfy two conditions. The first con- 
dition is 



• Ai » Aqcd- 

The matching factor / Wilson coefficient is calculated in perturbation theory. At a given order, 
the truncation error grows as /i decreases. Furthermore, if the operator couples to the Nambu- 
Goldtone (NG) boson due to spontaneous symmetry breaking (SSB), the effect cannot be corrected 
by perturbation theory and fails the matching. The condition is required to have these effects small. 
It should be noted, however, for some operators the SSB effect diverges as 1/w which has to be 
subtracted. Secondly, 

• /!< l/a. 

This is to reduce the discretization error, such like {pa)^. However, the {pa)^ > 1 region is used in 
typical calculations. It appears, in practice, jJ. <^ n/a is OK. 

These conditions provide the guideline to have the systematic error small for the RI/MOM 
schemes and the matching. A more detailed discussion follows by using an example of quark mass 
renormalization for the estimate of the systematic errors with rif = 2+1 DWFs. Here for simplicity 
we assume exact chiral symmetry. An issue of the small explicit chiral symmetry breaking will be 
discussed in Sec. 2A. The mass renormalization for m^ = Z,„m is conveniently calculated through 
the relation 

Z,„ = \/Zs = l/Zp\ (2.6) 



'Z„, = l/Zp 7^ i/Zs for Wilson type fermions. 
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Figure 1: Fractional difference (left panel) and average (right panel) of vector and axialvector vertex ampli- 
tudes as functions of [pa)^ for «/ = 2 + 1 DWFs at a^' ~ 1 .7 GeV [^]. Three masses are used to make an 
extrapolation to chiral limit. 



where Z„,, Z5, Zp are quark mass, flavor non-singlet scalar and pseudoscalar renormalization fac- 
tors. The ratio Zq/ZY is calculated from the traced projected vertex functions through the renor- 



malization condition Eq. ( |2.2| ), 



Ar = Tr[nr -^r]^^ 






(2.7) 



It immediately turns out that A5 7^ Ap due to SSB, i.e., Ap diverges as 1/m [Q |T^]. Therefore one 
needs to subtract the divergence or abandon the use of Ap. Here we proceed to take the latter choice 
and use ^ 



Ac 



1 

n 



Trps-/] 



Zs 



7 7 



(2.8) 



For Wilson type fermions, the pseudoscalar vertex must be used for the mass renormalization. The 
divergent pion pole contribution must be subtracted [|T2|]. Fitting the quark mass dependence with 
NLO OPE expression is one way. Or one can construct a object free from the pole combining the 
multi mass points [p^]. 

The scalar vertex, used in this example, does not have an infrared divergence. However, the 
naive unitary chiral extrapolation cannot remove a non-perturbative (NP) effect. If one has a handle 
on the direction of partially quenched mass, the leading NP effect can be isolated [^. Or one can 
study the response of the vertex function to the mass to evaluate its sensitivity to the low scale. The 
latter leads to an estimate of 7% systematic error for the present DWF case. 



The quark wave function renormalization is calculated from Ay through Eq. (2.3). With a 
chiral fermion formulation it is often averaged with the traced axial vector vertex A^ to gain statis- 
tics. In perturbation theory Ay = A^, but this can be violated at low momenta. The WT identity. 



i.e. the equivalence of Eqs. ( |2.3| ) and (|2.4|), holds non-perturbatively for the vector current. As 
the axial current couples with the NG boson, the equality is violated non-perturbatively. Fig. |] 
shows an example of the difference of A^ —Ay as a function of {pdf- . The difference gives an 



In quenched approximation there also be a divergence in the scalar vertex as l/m~ ||ll[|, which is suppressed in 
dynamical calculation from the quark determinant effect. 
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Figure 2: Left panel shows Z,^^^'-"^(jU) with ^^ — {paf', where cyan line indicates p^ = (2GeV)^. Right 
panel shows same data and those with RG running stripped off to get Z,^''^*^'^(2GeV) by LO to NNNLO 
PT as a function of matching scale. 



estimate of the NP contamination error. We will later use the window 1.3 < {pdf- < 2.5 where 
the largest difference is about 1 % which is small compared to the other systematic errors, such 
as the aforementioned 7%. The right panel shows the average (A^ + Av)/2. Now the quark mass 
renormaUzation is calculated through 



7RI/MOM 



(M) 



2As 



Aa+Av 



(2.9) 



where Za (= Zy) is obtained from a ratio of hadronic two point correlation functions of partially 
conserved and local axial vector current [^, Za = (j2^(x)f (0))/(Ao(.x)P(0)). 

The left panel of Fig. | shows Z^'^'^^{ii) with li'^=p^ as a function of {paf. Z^'^^ilGtV), 
for example, can be read off from the intercept with the cyan line drawn at p^ = (2GeV)^. This 
naive determination, however, might suffer from discretization errors, as well as the NP contami- 
nation error which has already been estimated. 

For the discretization error a common procedure is to strip off the renormalization-group (RG) 
running with perturbation theory. If the truncation error is negligible in a certain range of momenta, 
a linear extrapolation to {paY = from the range will remove a leading (pa)^ error. However, as 
a comparison between different orders of perturbation theory indicates, this extrapolation does not 
always give a better estimate than the simple method from the intercept (left). A particular problem 
of the extrapolation is that the remnant truncation error of the perturbative running gets enhanced 
by {pdf- — )• 0. Thus, if there is any indication of sizable PT truncation error, as is indicated in the 
figure, {paY — )• should not been performed. 

Fixing the matching scale pto & non-zero value could result in a biased estimate of the renor- 
malization factor due to a discretization error 0{p^a^) when only one lattice spacing is available. 
However, keeping the scale in physical unit for instance p = 2 GeV while taking the continuum 
Umit of the renormalized quantity will take away the leading {paY error automatically. This is 
clearly better than taking the {paY — )• limit at fixed a when a sizable truncation error of PT is 
suspected. 
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As the {paY dependence for Z,„ in the figure is flat with the NNNLO running, it does not matter 
whether to take the [paf — > or not in this particular case. The {paf — > value was adopted for 
2RI/MOM ^2GeV) in [^. Z^^(2GeV) is obtained by multiplying the conversion factor calculated to 



three loops [g|], Eq. ( [2.12| ). The truncation error is estimated from the size of the three-loop term, 
which turns out to be 6%. The resulting renormalization factor is Z,^^ (2GeV) = 1.656(157), where 
all the errors (mass dependence 7%, A^ —Ay 1%, PT truncation error 6%) have been summed up 
in quadrature with the statistical error. 

2.2 RI/MOM new developments 

As described above there are three competing effects for possible systematic uncertainty in the 
RI/MOM schemes which are: 

1. unwanted non-perturbative effects, 

2. {paY and higher, 

3 . a"+ ^ (/x ) and higher. 

The last error is relevant for the renormalization factor in MS, or in the case where the perturba- 
tive running are used for whatever purpose. There are several ideas to overcome one of them or 
aiming to reduce most of them. Here the methods are explained for the recent large scale calcu- 
lations for the quark mass renormalization by four collaborations: ETM for «/ = 2 twisted mass 



fermions with RF/MOM [|14|, |15|], QCDSF/UKQCD for «/ = 2 improved Wilson fermions with 
RF/MOM ^ [It]], JLQCD/TWQCD for rif = 2 and 2+ 1 overlap fermions with RL'MOM [|l|], 
and RBC/UKQCD for «/ = 2 + 1 DWFs with RI/SMOM and RI/SMOM^^ [|l^. 

For 1, ETM, QCDSF/UKQCD and JLQCD/TWQCD perform the subtraction. ETM and 
JLQCD/TWQCD exploit the partially quenched data for the fit to subtract the leading effects. 
JLQCD/TWQCD also performed the RI/MOM analysis on the £-regime where the SSB effect is 
suppressed. They demonstrated that A^ — Ay = holds well in the region of p^ where even a few 
percent difference was observed with p-regime data. RBC/UKQCD uses SMOM schemes (see 



Sec. 2.3) 



For 2, ETM performs 0{g^a^) subtraction for terms like {paY and T.^{pt/p^)ci^, obtained 
through analytic calculation in lattice perturbation theory [pO|]. QCDSF/UKQCD uses a similar 
subtraction, but, with a numerical integration in lattice perturbation theory. The two collaborations 
extrapolate remnant (pa)^ — )■ 0. RBC/UKQCD fixes p = 2 GeV, expecting the continuum extrap- 
olation of the renormalized quantity with a^ — )• removes the (pa)^ error. As a second method 
ETM performs the similar procedure as RBC/UKQCD and they checked the consistency with the 
first method after the continuum extrapolation. 

The RI/MOM to MS conversion factor for the quark mass is available to a,^ [|]. JLQCD/TWQCD 
estimates the size of the af contribution from the extrapolation of a" contribution for « < 3 and take 
it as the truncation error. RBC/UKQCD study the difference of Z^^ from two different RI/SMOM 
intermediate schemes for the truncation error of RI/SMOM — )■ MS perturbative conversion. They 
also check the consistency of the error with the size of the highest order in the conversion factor. 
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\q=PrP2 




Figure 3: A sample figure which makes 1/p^ contribution with SSB from Ref. [ph (left) and general momen- 
tum configuration of bilinear vertex (right) including the exceptional (q = 0) and symmetric non-exceptional 
(Pi — P2~ ?^) momenta from Ref. 



2.3 RI/SMOM schemes 



It was demonstrated in Ref. |Q] that the severe infrared sensitivity of the RI/MOM scheme can 
significantly be reduced by avoiding the use of exceptional momenta. The argument stems from 



Weinberg's theorem | |21| ] on the behavior of the vertex function for large external momenta, where 
a set of external momenta which has zero partial sum is called exceptional. 

In the case of bilinear operators with exceptional momenta, the momentum that flows in from 
one fermion leg can be rerouted to the other leg by one-gluon exchange (see Fig. |3|). If this happens 
the operator part allows a small momentum flow, which triggers SSB. This contamination only 
suppresses as l/p'^ at high momenta, hence could be sizable for the typical allowed momentum in 
simulations. This can be avoided by eliminating the exceptional channel by setting all p\, p2 and q 
non-zero in the right panel in Fig. ^. 

The RI/SMOM and RI/SMOM schemes are constructed for the quark bilinear operators 
in Ref. [p2|], exploiting the non-exceptional momenta and carefully ensuring the chiral Ward- 



Takahashi identities are satisfied. The symmetric configuration p^ 



-'^ — q^ is used, hence, "sym- 



metric" MOM scheme. A trial calculation [10] showed the RI/SMOM scheme is a promising 
alternative to the conventional RI/MOM scheme with reduced systematic uncertainties. 

The RI/SMOM scheme uses a different projector ^ oc (j[q^/q^ compared to ^ oc y^ in 



Eq. (|2.3|). The resulting quark wavefunction renormalization is equivalent to that of the RI7M0M 
scheme, i.e. z^^^^^^ = z^'™^'^. Another scheme named RI/SMOM^ uses the same projector 
as the RI/MOM scheme ^ oc y^ . It should be remembered, though, due to the use of the symmet- 



ric momentum, the resulting wave function renormalization is different, i.e. Zq 



,RI/SMOMv, 



/Z; 



RI/MOM 



Not only the reduction of the unwanted non-perturbative effect, but the better convergence of the 
PT for the matching to MS from SMOM schemes are observed. The matching was calculated to 



one-loop [|22|]. Evaluating the magnitude at each order in «« 



C„,(RI/SMOM^MS,At=2GeV,?i/ = 3) = 1 -0.015 + •• 
C™(RI/SMOM^^ -^ MS,;U = 2GeV,?i/ = 3) = 1 - 0.045 + • • 



(2.10) 
(2.11) 



0!.s (2GeV) = 0.2907 from four loop beta function has been used. 
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Figure 4: The difference Aj - Ap as a function of p^ GeV^ for a ~ 0.11 (24^) and 0.09 fm (32^) DWF 
lattices (left) [ p^ ] . Solid straight line is proportional 1 / p^ for guide for eyes. Right panel shows Z^^ (2GeV) 
with the SMOM or SMOMj,^ intermediate NPR scheme as function of [paf for the same a ~ 0. 1 1 fm DWF 
lattice. 



In the RI/MOM and RI'/MOM schemes the conversion factors are known to three-loop order 



C™(RI/MOM^MS,At =2GeV,«/ = 3) = 1 -0.123 -0.070-0.048 + ••• , (2.12) 
C^(Rl7MOM^MS,Ai=2GeV,«y=3) = 1 -0.123 -0.065 -0.044 + ■•• . (2.13) 

The convergence of the new SMOM schemes to MS is better. Moreover, it has recently been 



confirmed that this trend continues to two loops [24, 25], strongly suggesting the truncation error 
of the matching is significantly small for the SMOM schemes. 

The new SMOM schemes are being used for the quark mass renormalization in rif = 2 + \ 
DWFs with two lattice spacings by RBC/UKQCD collaborations ^^. The left panel of Fig. | 
shows the difference of pseudoscalar and scalar vertex amplitude Ap — A,, as an indicator of the 
unwanted non-perturbative effect, at the chiral limit. Not only is the difference finite, but small (1% 
at p = 2 GeV for instance), which is in clear contrast to the same quantity in the MOM scheme 
where it was divergent in the chiral limit as \/{mp^). Furthermore it exhibits the 1/ p^ behaviour 
expected from Weinberg's power counting and the possible symmetry breaking pattern for the non- 
exceptional momenta [^. 

The right panel shows the Z,^^(2GeV) for the a ~ 0.11 fm lattice, evaluated with the SMOM 



or SMOMy^ intermediate NPR scheme combined with one-loop matching [ ]22[ ] and two-loop MS 
running as function of the matching scale squared. The origin of the inconsistency could be the 
truncation error of the perturbation theory and the lattice discretization error. The non-perturbative 
effect is too small to explain this difference. The existence of the two independent estimates of Z^^ 
helps to estimate the systematic error. A preliminary estimate of the central value from an error 
weighted average at p = 2 GeV, taking into account the spread of the point around the linear fit, 
a variation with respect to the change of the matching scale is given as Z^^(2GeV) = 1.568(75), 
where the error includes an estimate of PT truncation error. This value is consistent with the same 
quantity obtained with the conventional RI/MOM scheme, Z^^(2GeV) = 1.656(157) [g], but with 
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a considerably reduced error. Further eiTor reduction is expected by incorporating the two-loop 
matching factors. 

The non-exceptional momentum configuration has been applied to the renormalization of non 
standard-model Bk operators, where the usual RI/MOM kinematics produces \/{mp^) type in- 
frared divergence. The non-exceptional momenta removed these divergences completely as ex- 
pected ^. 

An extension of SMOM kinematics to the standard-model four Fermi operator for Bk uses 
SMOM and SMOM^^ Zq. For the Y projector for the four quark operator, one can take the conven- 
tional type as well as the one with momentum transfer q^ in analogy to the SMOM and SMOMy^ 
projector. The 2x2 combinations make four SMOM schemes. Including the original RI/MOM 
scheme, five independent schemes lead to a better estimate of the central value of Z^^^(2GeV) and 
the systematic error [|I^]. 

The advantages and remaining issues of the SMOM schemes are summarized as follows. Some 
have already been mentioned. The RI/SMOM schemes improved from the RI/MOM scheme on: 

• Contamination of the unwanted non-perturbative SSB effect is greatly reduced. 

• Mass dependence is very weak. 

• PT matching to MS has far better convergence. 

• Signal/noise ratio for the non-exceptional momentum is improved from the exceptional. 

The second point is consistent with the first point and particularly convenient for «/ = 2 + 1 cal- 
culations, where one can safely keep the strange mass around the physical point. The last point 
is likely due to the insensitivity of non-exceptional vertex functions to the low energy fluctuation. 
There could be a room for an improvement on 

• enhanced effect of 0(4) breaking in momentum space. 

The last issue is visible in the right panel of Fig. § especially for the RI/SMOM scheme. It mani- 
fests itself as a non-smooth behaviour in (pa)^. Such an effect persists also in the MOM scheme. 



but smaller. The leading 0(4) breaking term appears at 0{g^a^) with a structure ofY,^pt/p^ [2C] 
for MOM scheme. As there is no precise knowledge of such terms in SMOM schemes, the effect 
has just been estimated from the spread of the data and taken into account in the systematic error. 
Given that now the two-loop conversion factor is available, this error could dominate the systematic 
error in the near future, then the solution to this issue will be desired. 

2.4 Other remarks in RI/MOM schemes 

For DWFs at finite fifth dimension Lv the axial Ward-Takahashi identity has an additional 
contribution to the naive extension from the continuum one from the pseudoscalar like the operator 
at mid point of fifth dimension. As such the partially conserved axial current £/" does not conserve 
in the chiral limit, which makes the renormalization factor of the current deviate from 1 [p^]. 
Denoting the additive renormalization of the quark mass in lattice units due to the non-vanishing 
midpoint term as m^es, the effect was estimated as Zj^ = 1 + O(mres) for the perturbative and Z^ = 



1 + 0(m^gj for the non-perturbative contribution |27]. The size of the correction was estimated to 

10 
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be no larger than 1% for ny = 2 + 1 DWFs at a ~ 0.1 1 fm. As the ratio of the local and partially 
conserved axial current is an input to the RI/MOM schemes to evaluate the quark wavefunction 
renormalization, this systematic error will propagate to the renormalization factors for multi-quark 
operators and the quark mass. This is not a dominant error in the current calculations. However, if 
the other errors are getting smaller, eventually this issue will need to be revisited. 

It is reported in Ref. [^ that the use of point sources for the RI/MOM scheme can result in 
underestimating the statistical error. The problem is enhanced towards the larger momenta. This is 
understandable because a larger momentum in the Fourier transformation of the sink position tends 
to sample a narrower region in the real space, thus effectively reduces the statistics. One solution 
is to have several distant source positions or to have a random source position. Another is the use 
of a volume momentum source, which effectively samples the entire volume. 

3. Developments in Schrodinger functional scheme 

Finite volume scheme with Schrodinger functional (SF) boundary conditions combined with 
a step scaling analysis provides a solid framework to perform a precision non-perturbative renor- 
malization [p8[]. The renormalization group running of the operators in the continuum theory can 
be calculated in a very efficient manner. Possible applications include the determination of the 
running coupling constant, renormalization of multi-quark operators, coefficients of Symanzik im- 
provement operators. The operator renormalization condition is applied to correlation functions of 
the operator in the bulk and distant boundary field(s). As such, the scheme is compatible for the 
Symanzik on-shell improvement of the operators. The practical difficulty is that the scheme needs 
tailored gauge ensemble generation for many parameter points for the step scaling and continuum 
limit. However, once the running of the particular operator is determined (with any gauge-fermion 
action), the remaining process is to perform the matching at given lattice spacing. 

The technique has been explained to some detail in the plenary presentations in the past con- 
ferences (see [^, ^]). In the plenary talk of Lattice 2002 by Sommer, who gave the last plenary 
presentation dedicated on the non-perturbative renormalization until this year, he listed three prob- 
lems which remain to be tackled. All related to the SF scheme: 

• Development of SF scheme for HQET at 1/m/,, 

• «/ = 3 simulation for non-perturbative estimate of as,v{pL), 

• applicability of SF scheme for the chiral fermions. 

By this year the answers of these questions have been given. The non-perturbative renormalization 
of HQET formulation at 1 /m/j is important to b-quark physics. In a recent lattice conference the 
calculation strategy has been explained in detail by Delia Morte [^. This year the Alpha collabo- 
ration reported on the calculation of ground state and first excited state spectrum of bottom-strange 
mesons, as well as the decay constant of B^, with completing the continuum extrapolation [p^]. 
Application to «/ = 2 improved Wilson fermion is underway. 

gO _ go mixing matrix elements have great importance to flavor physics (for example to con- 
strain the apex of the unitarity triangle of CKM matrix [3^, 31]). So far the renormalized operators 



are constructed in the static approximation to b-quark for «/ = 2 [ p^ ] and Uf = 2 [p3|]. If the 
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0{\/mh) operator is constructed, the systematic error from 0{\/m^) could be constrained to sub- 
percent level. However, it is yet to be realized. 

CP-PACS and PACS-CS collaborations have been performing «y = 2 + 1 simulations with 
non-perturbatively improved Wilson fermions. SF coupling ttsF is obtained by step scaling to give 
a high accuracy estimate of «/ = 3 coupling at a high scale, where matching to MS scheme is 
performed. Running down to charm threshold /i = nic by perturbative running to mach to rif = 4 



and then running up to «/ = 5, /i = mz, they obtained as{mz) [34]. See Sec. 4J for comparison 



with the other estimate. They push the project further to calculate the quark mass renormalization 



through SF scheme. The preliminary results are reported in [35]. The resulting renormalized quark 
mass is discussed in Sec. ^^ together with those from other discretization/NPR schemes. 

The SF boundary condition on the quark propagator naively interferes with the Ginsparg- 
Wilson (GW) relation at the boundary. As such, the GW relation has to be modified to properly 
define the SF formalism on overlap or domain-wall fermions (DWFs). The solution has been 



provided by Taniguchi for quenched overlap [|3g] and DW [ |37[ ] fermions through orbifolding. This 



method is successfully applied to the renormalization for the light quark mass and Bk with DWFs 



in the quenched approximation [g^]. It is extended by Sint [ |39| ] for the even number of flavors case. 



and Luscher [MOT for any number of flavors for overlap fermions. This year further development has 



been reported by Takeda [ ]41| ] on SF scheme construction in DWF with any number of fermions. 

The SF scheme is being used for twisted mass simulations as well. ETM collaboration tune 
the boundary counter term of "chirally rotated SF" for twisted mass fermions non-perturbatively 



[42] to ensure the automatic 0{a) improvement. The tuning has been performed for the quenched 



approximation. Application to mj = 2 is underway. 

4. Application specific developments/remarks 

In this section application specific developments in recent years and some comparison of dif- 
ferent schemes are provided. 

4.1 Strong coupling constant Us 

Due to the availability of several independent nf = 2 + \ simulations in recent years, lattice 
estimates of a^^{Mz) have now become very strong compared to the other methods. This has not 
been possible when only quench or «/ = 2 simulations were available, because of the inaccessibility 
of the threshold scale /i = m„ j or m, in perturbation theory. An estimate of as(/x) at rif = 3 for 
/I = 0(1) GeV allows one to use perturbative matching at charm threshold /x = m^ by perturbation 
theory (PT) to rif = 4, and then at /i = m^ to «j = 5, with RG evolution with anomalous dimension 
calculated with PT. The systematic error which arises in this calculation apart from the lattice non- 
perturbative part mainly comes from I) the PT machining of a particular scheme used in NP part 
to MS and 2) the ny = 3 o 4 matching at /x = m^ in MS scheme. Due to 2), one needs to select a 
scheme with which the matching is available to high order typically to a^. 

There have been four estimates of a^^{Mz) by the conference, 1) through ay calculated 
from Wilson loops [ p7| , p] ], 2) from moments of the charmonium current correlators [p3|], 3) from 



SF coupling [34], 4) from vacuum polarization function [^, 46]. For 1) one calculates Wilson 
loops to get ay {pi) through lattice perturbation theory, and then matched to a]^^(/x) at three-loop 
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Figure 5: Summary of a^'^{Mz) from rif ~2 + \ lattice simulations, compared with PDG 2008 average 
[p3]] (black). Results are from Wilson loop [E3] and charm correlation functions [45| with staggered (red), 
SF coupling [B4| with Wilson (green), and vacuum polarization function [M with overlap fermions (blue). 



(NNNLO). 2) uses the continuum three-loop expression (NNNLO) of a moment of charm current- 
current correlation function. This is a by-product of their charm quark mass estimate for which 
they use same correlation functions. 3) use step scaling of the SF coupling and matching to MS. 
Although this has been obtained only to two loops (NNLO), the matching is performed at high scale 
IX ^ trie, where the PT truncation error is small compared to the other errors. The renormalization 
scale is set from the inverse linear lattice extension jx = \/L. Once tts (m) has been obtained, it 
is run down to /x' = lUc to match to nf = 4. This matching and running are done with NNNLO. 4) 
the continuum vacuum polarization function has been obtained through operator product expansion 
and the relevant coefficients has been calculated to NNNLO. The renormalization scale is set from 
the size of the injected momentum at the current. A compiled summary of a^^(Mz) estimates is 
shown in Fig. |5[ Results from different method / lattice discretization are consistent with each other 



and with the PDG2008 average [yj. 

All of these calculations reached a remarkable precision of the strong coupling constant. One 
has to keep in mind, however, all of them uses the perturbation theory at /i ~ m^, where non- 
perturbative effect might still be important. Analysis on ny = 2 + 1 + 1 simulations on finer lattices 
will be able to address this issue. 

There are several technical development on determining Aj^ using the Landau-gauge ghost 
and gluon propagators, which has been applied to «y = and 2 [^, 50, 51]. 



4.2 Light quark mass 

A comparison of the light quark mass from different lattice discretization and renormalization 
schemes is presented by Scholz [ |60| ] in this conference. Here for the purpose of discussing the 
renormalization issues the strange mass compilation is made with the available data in the literature 
(Fig. i). 

For rif = 2 + \ dynamical fermion computations, a remarkable precision has been achieved 
by HPQCD collaboration [^, 52] who exploits a precise estimate of MS charm quark mass from 
the charmonium current correlator and the bare mass ratio ms/m^. The use of highly improved 
staggered quark (HISQ) and the fine lattice spacing made it possible to use the same action for 
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Figure 6: Summary of strange mass m^^(2GeV) for «/ = 2 and 2 - 
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both light and charm quarks. This method does not require the conventional lattice renormaliza- 



tion procedure. RBC and UKQCD collaborations [53] have used RI/SMOM schemes described in 



Sec. 2.3 on DWFs with two lattice spacings for the continuum extrapolation. The MILC collabora- 



tion [ p4| ] uses Asqtad action with two loop lattice perturbation theory [ |62| ] for the renormalization. 
These three results in the continuum limit agree with each other. One cannot find any issue for the 
perturbative vs. non-perturbative renormalization here. 

For results from single lattice spacing, Lytle reported the strange mass with MILC coarse 



lattice with RF/MOM NPR []55[]. He compared with the strange mass with perturbative renor- 
malization on the same lattice and found NPR gives about 20% larger value. PACS-CS reported 
the calculation from improved Wilson fermions directly simulated on the physical ud and s point 
with Schrodinger functional NPR [^,57], which has become 30% larger from their earlier estimate 
with perturbative renormalization referred in Scholz's compilation []60|]. These differences between 
the perturbative vs. non-perturbative renormalization are not necessarily inconsistent with the ob- 
served consistency for the values after the continuum extrapolation. This is because the difference 
could arise from the discretization error in different renormalization schemes. JLQCD/TWQCD 
has the estimate from overlap fermions with RI/MOM NPR. For these three studies, extension to 
the second lattice spacing are anticipated. 

There has been several Uf = 2 calculations reported for the light quark mass. In this figure 



of Ms compilation only earlier studies with Wilson fermions by Alpha [ ]58[ ] and QCDSF-UKQCD 



[ p9[ ] collaboration are shown. The new calculation of average ud quark mass in the continuum limit 



with twisted mass fermions combined with the RI'/MOM scheme has been reported [|63|, |64p. The 
result is consistent with rif = 2 + 1 continuum values 



4.3 Three quark operators 

A class of three quark operators which are a part of baryon-number violating four-fermion 
operators induced by OPE from grand unified theories are important for the estimate of proton 
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lifetime. The operator classification, construction of RI/MOM scheme, and one-loop matching to 



MS have been carried out in Ref. [65|. 



The RI/MOM scheme of more general three-quark operators are constructed by QCDSF- 
UKQCD 1 66], for the proton decay matrix elements and nucleon distribution amplitude [^]. 



4.4 Static heavy quark 

As mentioned, the SF scheme has been developed for b-quark system namely heavy-light bi- 
linear and AF = 2 four quark operators with the static approximation. Application of the RI/MOM 
scheme to this system was tried long ago [l68|]. The subtlety to have a mass independent scheme 
for the static quark formulation, where the self-energy has power divergence, prevented further 
development. In this conference, Izubuchi showed his study on the static quark self-energy, hoping 
to resolve this issue by making use of the odd part in energy of the Fourier transformed static quark 
propagator. Further investigation/demonstration is desired. 

4.5 Relativistic heavy quark 



The Fermilab-Tsukuba approach [ ]69| , [70| ] the heavy quark uses relativistic clover type fermions. 
The (m/ja)" error with all order in n is removed using perturbation theory. The Columbia group re- 
fined the strategy of this approach to correctly count the independent operator coefficients [ JTll ] and 
developed the non-perturbative determination with step-scaling [^]. Later another tuning proce- 
dure was proposed with by-passing the step scaling and matching experimental spectrum directly. 



It has been applied to charm []73[] and bottom [72] systems 



5. Related development in lattice perturbation theory 

Lattice perturbation theory has been developed for some quantities. For the quark mass renor- 



malization, two-loop results are available for staggered [ ]62[ ] and Wilson(W)/clover(C) fermions 



[74]. The latter calculates the renormalization constant in RI7M0M scheme. Once that is obtained 



the continuum matching can be used to convert it, for example, to MS, where the conversion factor 



has already been calculated up to three loops [|8| [23[ ]. 

0{g^a^) calculations have been performed on the quark self-energy and bilinear operators 



for W/C [EQ], twist two bilinears for W/C/twisted mass (tm) [75], four Fermi operators for W/tm 



fermions [76]. The first one has been used in combination with the non-perturbative RI7M0M 



renormalization to have 0{g^a^) effect on p^ and Y^uPt/p^ subtracted (see Sec. 2.2). 



Numerical stochastic perturbation theory (NSPS) pl\ \ is being developed for the improved 
gauge and Wilson fermion action on the quark self-energy [ |78| ] at three loop through the RI7M0M 
renormalization condition imposed in the lattice perturbation theory. Two-loop result for the 
Symanzik improvement coefficient csw and ca for Wilson fermions has been obtained [[7^]. The 



Landau-gauge ghost propagator has been calculated to three loop [|80|], anticipating the application 
to estimate of a^. 

6. Concluding remarks 

Today, the nf = 2+ I simulations, incorporating the full dynamics of three lighter quarks 
that cannot be accessed from the perturbation theory, are available with four different fermion dis- 
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cretizations. Thus various crosschecks on the calculated quantities can be performed. Furthermore, 
Hf = 2 + 1 + \ simulations, incorporating the dynamical charm effect, for two fermion actions are 



being carried out [|81|, |82p. In these circumstances the renormalization in lattice QCD has become 
more important than ever before. Non-perturbative renormalization, which was reviewed here is a 
powerful tool to tackle the precision calculation of the fundamental parameters of standard model 
and hadron-operator matrix elements. In recent years, the Schrodinger functional scheme has been 
applied to the strong coupling constant and light quark masses for nf = 3, non-perturbative tuning 
of the parameters for the simulations, as well as the b-quark systems. Due to the window problem, 
RI/MOM schemes at the typical lattice spacings of current simulations could suffer from larger 
systematic errors than the SF scheme. Yet, there has been several ideas of bringing the error under 
control, which have been discussed in this review. As a promising new direction the RI/SMOM 
schemes were described in some detail. A small summary was made for the development on the 
lattice perturbation theory, as well. The comparison of the strong coupling constant and the strange 
quark mass estimate in the continuum limit have shown consistency despite of the difference of the 
lattice action and renormalization method used. The success stands upon the cultivation of not only 
one but multiple relevant technologies in the lattice QCD community over the last two decades, and 
would motivate further improvement/development in this stimulating field. 
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